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ABSTRACT 

The use of an electronic analog computer in solving the differentisl 
equation of motion of the free piston gas generator has been investigated. 
To accomplish this end, the free piston has been treated as a specialized 
mass~spring system, wherein the forces were derived from thermodynemic 
considerations rather then from elastic considerations. The first law of 

/ 
thermodynamics was solved on the computer three concurrent times in order 
to produce electronic analogs of an internal combustion engine, a recipro- 
cating compressor, and a gas spring. Coulomb end viscous friction were 
also simulsted. The results of these four analogs were then combined and 
integrated as indicated by Newton's second law. The computer yielded a 
solution which represented the piston motion, A total of 26 amplifiers, 
seven diode limiters, and six electronic multipliers were used to solve 
the problen, 

The characteristics of an existing free piston machine were intro- 
duced into the computer as a final check of the analog. The accuracy of 
some of the results was less than desired due to errors in a two-quadrant 
electronic multiplication involved in the gas spring analog. The other 
three analogs performed very setisfactorily. For accurate and repro- 
ducible results, some other method, such as direct function generation 
or single quadrant multiplication, must be used in the simulation of 


the gas spring. 
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1. Introduction, 

A great deal of effort has been put into the solution of the free 
piston engine problem, All of the methods currently in use are digital 
methods involving a point-by-point calculation of the variables asso-= 
clated with such a device, These methods are time-consuming and cumber- 
some due to both the trial and error element and the many geometric and 
operating parameters involved. In addition, the digital methods can be 
made tractable only by incorporating many numerical assumptions; some 
of which introduce considerable error. 


cA) has developed a "Thermodynamic-dynamic analysis" whereby 


London 
the calculated parameters of a "standard" machine may be scaled by 
certain affinity relationships to establish the design of larger or 
smaller machines, This method eases somewhat the numerical burden of 
calculation, but the results proved to be quite sensitive to the assump= 
tions made; and in order to predict the characteristics of an existing 
machine, the assumptions were relaxed in four successive groups before 
the desired accuracy was obtained. This procedure still presented a 
formidable numerical task, 

In addition, a solution by digital methods becomes more complex in 
the investigation of part-load characteristics, 

The purpose of this investigation is to determine the possibility 
of solving the differential equation of motion of the free piston on 
the Boeing Electronic Analog Computer (BEAC), which is described in 
Appendix VI. If this can be done, the computer will electronically 
Simulate the free piston, presenting as a solution the piston motion, 
with most of the pertinent quantities also available. The advantages 
of this method include: 


1, The feedback principle, typical of analog computation, 
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eliminates trial and error, 

2e <All of the parameters may be taken into consideration, and 
the effect of each may be observed. 

5. The solution is immediate, thus permitting the study of many 
combinations of parameters in a short time. 

4. Fewer simplifying assumptions are required, and SAE may be 
relaxed to any degree compatible with the size of the computer, 

5e Parteload behavior can be studied as easily as full-load 
behavior, 


6. Time may be scaled to permit slow=-motion studies, 





2e General Background of Free Piston Systems. 

Figure 1 illustrates schematically two of the possible free piston 
configurations. Figure l(a) illustrates the free piston air compressor, 
where the useful output is compressed air; and Fig. 1(b) illustrates the 
free piston gas generator, where the useful output is high pressure- 
temperature gas for use in a turbine, 

This investigation will encompass only the free piston gas generator. 
A slightly larger computer is required to simulate the free piston com- 
pressor. 

The free piston machine depicted in Fig. 1(b) consists of two 
opposed pistons, the inboard ends of which operate as a two stroke diesel 
engine. Another part of each piston acts as a reciprocating air com= 
pressor, As the pistons move outward, the expansion work of the ane 
and compressor is stored as energy in the bounce chambers, This energy 
is then used on the inward stroke for the compression and delivery of 
the air in the compressor and the compression of the engine charge. 

The outstanding features of this free piston machine are: 

1, Simplicity of construction due to absence of orank 
mechanisms, 

2e Opposed piston configuration eliminates vibration. 

Se No piston side thrust arising from crank inertia forces, 

In the case of the gas generator, it would be well to review the 
features of a typical cycle, 

Referring to Figs 1(b), which illustrates, say, the pistons just 
starting the instroke. Beginning at this point, the following successive 
processes may be visualized: 

1, High pressure in the bounce space forcing the pistons 
together; compressing and delivering the air in the compressor space, and 
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Figure 1. Typical Free Piston Configurations. 





compressing the engine oharge when the scavenge ports are oovered. 

2 Reversal of piston motion occurring during injection and 
burning of fuel. 

3. The expansion of the fluid in the engine cylinder and the 
intake stroke of the compressor; while the bounce cylinder is storing the 
expansion energy of the engine and compressor. 

4, Opening of the engine ports during the latter portion of 
the outward stroke, delivsring working fluid to turbine, 

5. Pressure in bounce space causes reversal of piston motion, 
and cycle repeats. 

The preceding cycle is represented on Temperature=Entropy coordi- 
nates in Fig. 2. It is to be noted that the processes involve different 
masses, and can be represented only on a "per pound" basis, This diagram 
is similar to that in common use in free piston literature. It should be 
pointed out that the following two assumptions are made in producing this 
diagrams 

1. The scavenge air is furnished at constant pressure and 
temperature, related to the Seer gees discharge state, 

2e Delivery to the turbine occurs at constant pressure and 
temperature, 

These two assumptions are intimately connected with the mechanical design 
of the engine ports and the physical size of the turbine receiver. It 

is beyond the scope of this investigation to oonsider such problems. The 
foregoing, therefore, constitute the first two assumptions in this 
investigation. It must be pointed out, however, that these assumptions 
must be relaxed if it is desired to investigate transient behavior, inas- 


much as these particular parameters are the primary transients. Before 
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Process 1-2 - Compression in Compressor, 
Process 2-3 - Pressure Drop into Fngine,. 
Process 3-/, - Compression in Engine, 

Process 4-5 = Burning of Fuel. 
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Process 6=7 - Blowdown as Exhaust Ports Open. 
Process 3,7=-& - Mixing and Scevenging. 
Process 8=-9 - Expansion in Turbine, 


Figure 2. Temperature - Entropy Representation of Free Pistcm 


Generator-Turbine Systen. 
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proceeding further, the following are some additional assumptions” 
covering the device as a whole: 

1. That the compressor intake is at 70°F. and 14,7 psia. 

2. That the operating fluid obeys the ideal gas laws, 

Se A nominal pressure drop of five per cent across each valve 
and engine port. This is arbitrary, and may be refined to any desired 
value, 

4, That the geometry is symmetrical with respect to a plane 
passing through the center of the machine perpendicular to the hori- 
zontal axis. The end projection of this plane will be referred to as 


the "center line," 


*Refer to Appendix I for all of the assumptions used. 





3. The General Analog of the Free Piston. 

From a dynamic point of view, the piston in a free piston machine may 
be considered as a rather special sort of mass-spring arrangement, wherein 
the "springs" yield forces derived from thermodynamic considerations 


rather than from elastic considerations. 


Li if 
NN I / 
SAA 


| 


Figure 3. Spring Analogy of the Free Piston. 


The system may be visualized as similar to that depicted in Fig. 3 above, 
where the spring forces are proportional to the conventional Pressure- 
Volume representation of the appropriate space. Not shown, but also 
present, is the force of static and viscous friction. Qualitative 
representations of these forces appear in Fig. 4. 

To illustrate the analog computer approach to this type of problen, 
consider the simple mass-spring arrangement illustrated in Fig. 5, The 


differential equation of this motion may be immediately written down as 
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where m ig the mass 
v is the damping constant 
k is the spring constant 
7 is the acceleration 
J is the velocity 
fis the displacement 


By solving for the highest 
derivative, 


4--(K7+x ?). 





Figure 5. Simple Spring-Mass System. 


Recalling that the operational amplifier can add and/or integrate 
with respect to time, a computer schematic may be immediately devised to 


solve this problem. This is illustrated below in Fig. 6.” 


y [> 


B\k< 





Figure 6. Electronic Analog of the Simple Spring-Mass System. 


"Refer to Appendix II for analog computer table of symbols, 


10 





This simple problem points out the utility of analog feedback. That is 
to say, the immediate stabilization to the correct solution by using this 
solution to generate the derivatives, which are then integrated to form 
the solution again. 

Turning now to the free piston problem, and referring to Fig. 5; the 
equation for this motion may be formulated in accordance with Newton's 
second law as 

MX = FRX) + Fo(%,X) - Fo(%,X) - Pe( Zz), ou 
where M is the piston mass 
Fg ‘7,%,) is the force on the engins piston 
Fo(z,X) is the force on the compressor piston 
F,(4,%) is the force on the bounce piston 
Fe(%) is the force due to friction 
is the piston displacement 


is the piston velocity 
is the piston acceleration 


asd. Ps 


In a manner similar to the simple mass-spring problem, a computing ar- 
rangement can be immediately proposed. Suoh an eer ee outlined 
in Fig. 7. In this figure the block representing the engine, compressor, 
bounce, and friction, are rather complex computing networks with many 
more inputs than shown, The overall principle, however, is basically 
the same as that outlined in Fig. 6. 

This fact has two very important consequences in that : 

1, A completely compatible set of parameters must exist in 
order for the computer to produce a sustained oscillation, If this cone 
dition is not met, a situation exists which is analagous to overdamping 
or negative damping of the simple system, 

2e Sinoe equation 3,1 also expresses the work-eenergy equality, 
a solution of this equation will automatically satisfy the condition 
that the engine work appears as the sim of peer assor’ and friction work, 


The computer arrangement as shown in Fig, 7 will produce displace- 


ll 
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Figure 7. Basic Computer Arrangement for the Free Piston 
Problem. 
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ment as a function of time, Also available will be the various foroes 
or pressures contributing to the motion. Related items such as air-fuel 
ratio, scavenge ratio, and powor output will not be available direotly 
from the computer but oan be easily caloulated, 

The inputs to the system will be : 

1, Geometric parameters, which will be eiked for a given 
machine, 
2e Initial pressure conditions, which are: 
& Intake and delivery pressure of the nencuanteer « 
be Pressure at the start of engine compression. 
ec, An initial pressure in the bounce BPACGs 
Se Amount and heating value of fuel, 
4. Thermodynamio parameters such as polytropic exponents, 
and ratio of specific heats, 

All of the above parameters may be varied to suit almost any oon- 
dition compatible with the assumed hawt, For instance, part load 
behavior may be simulated by merely changing the appropriate pressure 
initial conditions, and then determining the range of fuel energies for 
which the solution is stable. Of course, a different amowumt of fuel 
applied to the same pressure conditions will lead to a different dis-~ 
placement-frequency solution, 

The usefulness of the analog computer approach to this problem is 
that: 

1. The parameters are introduced as various constants and 
initial conditions and are not manipulated at all, 

2e The computer solves the differential equation of motion. 

5. Other quantities may be determined by a minimum amount 
of caloulation,. 
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It is to be re-emphasized that the analog computer eliminates the 
trialeand-error elements from this problem if the geometry is used as a 


starting point. 
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4. Geometrio Parameters, 
There are eight primary geometric parameters and variables asso0- 

ciated with the analog approach to the free piston problem. These ares 

1, Area of engine piston, Ap. 

2. Area of compressor piston, Age 

3. Area of bounce piston, Ap. 

4. Engine piston displacement from centerline, x. 

5. Compressor piston displacement from head, y. 

6. Bounoe piston displacement from head, z. 

7. Engine piston position at port closure, Xpe 

8, Piston mass, M. 


Referring to Fig. 1(b), it may be readily seen that 


Ay a Ap + Age 4.1 
y s x + Cy 402 
and 


where c and b are the respective clearances of the compressor and bounce 
cylinders when the engine pistons are touching. Since the piston is in 


one piece, it is obvious that 
dy. @3 4.4 


These relationships are necessary in that x, y, and z and their time 
derivatives are used in the generation of the appropriate forces, 

The mass, M, may be taken as any reasonable value compatible with 
piston size. A computer solution indicating a piston density greater 


than that of steel must be rejected. 
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The treatment of the free piston generator will encompass one-half 
of the system, which will be the right hand piston as illustrated in 
Fig. 5. Positive forces will be considered to be acting to the right. 
Since one of the computer processes will be the summation of these forces, 
the associated sign inversion requires that the negatives of the engine 
and compressor forces be generated. Since one-half of the system is 


being treated, all extensive quantities must be multiplied by two, 
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5. Thermodynamic Parameters. 

Inasmuch as the computer is solving a mechanical problem, all thermo- 
dynamic quantities will be assumed to be given in compatable mechanical 
units. The use of Joule's constant will always be implied whenever 
necessary. Energy leaving a system as heat will be considered negative. 
5.1 Pressures. 

The pressure taken as the fundamental variable will be the turbine 
receiver pressure, P;. This will be assumed to remain cyclically constent. 
Using this as a starting point, it is possible to work backwards through 


the system to establish all of the other pressures. These are 





ie 

Fe re 2 5 lel 

P= te . Fe 5.1.2 
rs keh ? 

esees 3 ee 5203 
a rels ra 

and 
P. «x lol 5.1.4 


where Pp, is the engine pressure during scavenging, 
Pz is the pressure in the scavenge cylinder, 
P, is the compressor delivery pressure, 
Ps; is the compressor intake pressure, 
P, is the ambient pressure, and 
r is the ratio of downstream to upstream pressure for 


"e" —~ Engine exhaust ports, 

"ag" ~ Engine scavenge ports, 

"qd" —~ Compressor discharge valves, and 
"i" ~ Compressor inteke valves. 


5.2 Volumes. / 
In all cylinders, the volume is the product of total displacement 

and area, 

5.3 Temperatures. 


\ 


Temperature will not be used as such, but rather indirectly in the 
= / 


determination of the mass of air charged to the engine and the power 


output of the machine, 
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5.4 Mass of Engine Charge. 

The mass of air charged to the engine is not at all a necessary para~ 
meter insofar as the computer is concerned, but must be known in order to 
determine the air-fuel ratio. 

Using the assumption that there is no heat gained or lost in the 
scavenging space, it can be postulated that, upon completion of the sca- 
venging process, the engine charge temperature, Tp» will be essentially 
that temperature of the compressor discharge. The mass of engine charge, 
Mo, Will then be 

Me = Po Vp/RP ; 5e4o1 
where Vp is the engine volume at the point of port closure. 
Applying the compression ratio of the compressor to the above equation 


yields 
_PeVp 


where T, is the ambient temperature 
n, is the polytropic exponent of compression in the compressor. 
This result may be expressed in terms of Py by applying relationships 
Sell and 5.1.6: 
( nel ) 


(a-Ne) 
ian Vp fe beg ah L. 6.403 


Ro te 


If all r's are assumed the sam, the above equation reduces to 


-Te) _ (n,-4 : 
Me = at a ees 





52404 


¢ 


Assuming that r is 0,95, and further postulating that Ny is 1.50, the 
above relationship becomes 


0.7 
i ; 5.4.$ 


mefy = 0,0113 P 
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where Py, is expressed in pounds per square inch. For convenience, 
equation 5.45 is presented graphically in Fig. 8. 
5.5 Power Output of the System. 
To determine the power delivered to the turbine, an overall energy 
balance is necessary: 
mah, + me HV = (mp + mg)he + Qrejs DeSel 
where mq is the mass of air delivered by the compressor per cycle. 
mp is the mass of fuel injected per cycle. 
h, is the enthalpy of entering air. 
ht is the enthalpy of exit gases. 
HV is the lower heating value of fuel ; : 
Qrej is the energy as.heat leaving the system via the redding 
water per cycle. 
Equation 5.5.1 makes the customary assumption that the sensible enthalpy 
of the fuel may be neglected, 
The quantity mg may be readily determined by noting from the come 
puter solution the length of the intake portion of the compressor stroke, 


and applying the ideal gas law 


7 


= PyAgY¥i/ RT, * 6.5.2 
where Y; is the length of the intake portion of the compressor stroke. 
The quantity mp HV is a known operating variable and hp, is tabulated. 
The quantity Qre4 is composed of: 
1, Energy as heat leaving the engine. 
2e Energy as heat leaving the compressor. 


Se nergy lost through frictional heating. It will ‘be assumed 


that this will be totally lost to the coolant, 


The determination of Qrej Will be deferred until more specific 


differential equations are formulated, c ‘ 


Assuming for the moment that Qrej is known, the only remaining une 
known in equation 5.5.1 is ht, the enthalpy of the exit gasses. From 
\ 
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Figure 8, Engine Cherge Mass per Cubie Foot Pistom Displacenent, 
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this quantity, the temperature of the exit gas may be determined, and 
also the ideal isentropic work available per pound of fluid. The pro- 
duct of this work, the total mass delivered per cycle, and the frequency 
will then give the air horsepower of the system. It is to be noted again 
that this treatment involves one-half of the machine, so that the "total" 


mass mentioned above is twice that given in equation 5.5.2. 


2) 





6, Electronic Analog of an Internal Combustion Engine. 
6.1 Generel. 

The force exerted upon the engine piston may be most appropriately 
represented by referring to a typical indicator diagram such as depicted 
in Fig. 4(a) on page 9. Since the free piston device has no crank mech- 
anism, the only fixed geometrical point is that of port closure. The 
thermodynamic stete at this point will be kmown if it is assumed that 
successive cycles are identical. All other point of force as a function 
of displacement are determined by the dynamics of the system. In order 
to simulate this situation it is necessary to obtain a computer solution 
of the differential equation of the first law of thermodynemics, 

In the following analysis, the approximate air-standard is used, which 
assumes 2 

1, Air is the working fluid, and 
2. The mass of fluid within the cylinder remains constent. 
The latter assumption is also derived from the following two 
practical considerations: 
1. The mathematical complication involved in a changing 
mass system, while possible to solve on the computer, would require 
more components than are available. 
2. The mass change is within the same order of magnitude 
as the accuracy of the computer, this there would be no observable 
gain from such e refinement. 
6.2 The General Differential Equation for thc Internal Combustion Engine, 
The first law of thermodynamics may be written in differential 


form as 


AQ = ol W ola 5 6 .omn 
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where dQ is the net energy entering or leaving the system as heat, 
aW is the work done by the system. 
dE is the change in internal energy of the system. 
By taking the time derivative of this equation, and by replacing dW and 
dE by their respective equivalents in terms of pressure and volume, the 
following equation is obtained. 
AQ AV | dv ap) 
——= = a + 
at ~ Par a K-42 Pat Vat . 66c0on 


By solving for the highest derivative of P35 
dP a [yy -4)92 - w Pav 
oe fie) K Poe . 


dt ~~ V 6.2.5 
Now, the volume is given by 
V = Az Xx J Belted 
and 
V lie 3 
dt Ed Ag Bib 


When equations 6.2.4 and 6.2.5 are substituted into equation 
6e2e2, the final result is 
Ale al | Sie eae KP. x 
dt XA] Ace dt 
where Pp is the engine pressure, 
This equation may be solved on the BEAC in the following steps. 
1. Integrate dPr to form -Pp with the appropriate initial 
dt 
condition, Poe 


2, Multiply -Pg by kx to form -kPgx. 


: K-l dQ _ | 
5. Add k=l dq to form [qt 48 KP. 2}. 


AR 
4. Divide this sum by x. This result is, by equation 6.2.6, 


dPr, which is then fed back to Step l. 
at 
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The above steps are represented schematioaldge in Fig. 9. This dia- 
gram also shows considerable computer detail which is explained in the 
following notes. 

1. The inputs to Amplifier 5 are to separate the oan into 
"fuel" and "non-fuel" terms in order that these quantities may be varied 
independently. In addition, the "fuel" function is non-zero for only a 
portion of the cycle. The output of this amplifier gee net. 

2. Tho relay in the feedback of Amplifier 1 establishes the 
initial condition, Po» whenever x is greater than Xpe 

Se The oombination of the electronic multiplier and Ampli- 
fier 2 is pre-set to give one-fiftieth of the produot. 

4. The multiplier in the feedback of Amplifier 4 has the 
effect of dividing by x, where-x must always be negative for stability. 
The output of this amplifier is pre-set to give ten times the quotient. 

Se Amplifier 3 is set to give ten times the sum of - a(t A %) 
and 70 (Ge ey in order to compensate for the constants inherent in 
the multiplication and division. 

To gain a better insight into the computer solution of equation 
6.2.86 a few simple qualitative examples will be disoussed in the followe 
ing three sections. 

6.5 Computer Simulation of a Reversible Adiabatic Process, 
The reversible adiabatic process is characterized by equal to 


zero, in which case equation 6.2.6 reduces to 


d Fe => = ares 
dt x . Oe et 


It may be easily verified that Pee constant is the solution to the 
above equation, which defines the familiar isentropic process. The 
oomputer solution to equation 6.3.1 appears in Fig. 10 with various 
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Figure 10, Computer Simulation of an Isentropic Process, 





initial conditions, The input "x" was a displaced sinusoid simulating a 
piston motion at about one cycle per second. This was produced from the 
arrangement shown in Fig. 11, where 2 is an auxiliary variable solving 
the equation Zz + ur3 =O, Also shown are the "y" and "z" variables 
which will be used in the forthcoming discussions of the compressor and 
bounce cylinders, This is an artificial stroke input and will be by= 
passed when the entire free piston analog is assembled. sa 

The solutions depicted in Fig. 10 can be seen to be excellent for 
all but the lowest values of initial condition. The plotted points are 
the numerical values of the indicated equations, and the solid lines are 
the computer solutions. In all cases, the integration was stopped when= 
ever "x" was greater than 30 volts, The departure from the correct 
solution for the low value of initial condition is due to the fact that 
the inherent electronic errors are magnified whenever the computer is 
working with small voltages, This condition must be avoided if at all 
possible, 

It is to be pointed out that the multiplication of =Seey x 
involves two quadrants, That is to say, =Pp is always negative, while 
x is both positive and negative, Care must be taken to insure that 
this multiplication yields the same result in both quadrants. A slight 
maladjustment in this multiplier will not return the solution to the 
starting point since dPy in one direction will be different from dP, in 
the other, a kal 


6.4 Computer Simulation of a Reversible Polytropic Process, 


The reversible polytropic process may be achieved by letting 


<e = “fe (x-» Ax, 644.1 
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Upon substitution, equation 6.2.6 becomes 


ae _ _ nlex | my 
At x . e@*@ 


It may be easily verified that Pex” = constant is the solution to this 
equation. It should be pointed out that while eS ae is as 
defined in equation 6.4.1, the same result may be achieved on the computer 
by letting dQ be zero, and applying the coefficient n instead of k after 
Anplifier ee suggested by equation 664e2.6 ° : 

Figure 12 illustrates the computer solution of a few samples of this 
case, again with excellent agreement. At this point, it is important to 
realize that there is no net loss of energy from the system over a cycle. 


The energy lost during compression is restored in the same manner during 


the expansion, or 


(Se) 48 = 0 6.4.3 


This may be easily verified on the computer by integrating equation 

604.1. Figure 15 shows displacement, pressure, dQ , and \/dQ\y, as 
dt dt 

functions of time. .As expected, Hague zero at the end of 4 cycle, 


~ 
6.5 Computer Simulation of a Nonepolytropic Process. 


Consider a process which is characterized by a net loss (or gain) 
of energy. This may be accomplished by touting 2G some funotion other 
than that defined by equation 6.4.1. ‘Suppose, for example, that dQ 
were a small negative constant. The resulting differential is 
not easy to solve, However, it can be deduced that a cyclic change in 
volume will not return the pressure to the starting point, since aPE 
is different for each direction of motion, A few computer eLueiels of 


processes of this type are shown in Fig. 14. With the exception of the 
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isentropic case, there is net negative work done. This is, of course, the 
work equivalent of the energy lost by the system as heat, Figure 16 
illustrates the same solutions as a fumction of time. Here it is noted 
that the net energy, Q, appears as a ramp function. The discontinuity in 
the pressure curve is due to the operation of the initial condition relay 
when the displacement exceeded 30 volts, 

This and the preceding examples serve no purpose in the overall aim 
of this investigation, but are presented in order that the reader may 
become familiar with the rather peculiar sort of mathematical manipula- 
tions used in setting up the analog computer. Im addition, they may serve 
as useful check points in the event that this problem is set up on a 


different computer, 
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7. Electronic Analog of the Standard Diesel Cyole. 
7,1 General 
Equation 6.2.6, and the associated computer arrangement illustrated 
in Fig. 9 represents the most general case of an engine cycle, wherein 
any form of heat removal or addition may be simulated by the appropriate 
selection of functions for the dQ terms. The investigation of solutions 
of this type, encompassing various methods of heat removal and fuel 
injection, is a very broad subject, and should be treated as a separate 
investigation. It will not be practical to do this within the wider scope 
of the free piston problem. In order that the results from this investi- 
gation may be more readily analysed, a standard diesel cycle, will be 
employed for the engine, This has the additional practical advantage 
that excessive computer voltages will be avoided, 
7.2 The Standard Diesel Cycle, 
The standard diesel cycle consists of the following processes3 
1. Polytropic compression, at exponent n, of the ongine charge. 
Ze Constant pressure addition of energy over a portion of the 
outs troke, 
3. Polytropio expansion of the engine charge and combustion 
products, at exponent n!', 

4. Constant volume heat rejection to exhaust, 

5. Constant pressure scavenging process, 
Referring to equations 6.2.6 and 6.4.1, these processes may be desoribed 
mathematically as follows$ 


1. Compression process. 


Ale = = alex Tet " 
t >) eGe 
d n ro 
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and Lo 
G2) = Ac (tae) Re. ~ 


2. Combustion process. 


Af 
(24). ie Tebed 


and 


Se Hxpansion process, 


_ Ged ae fe % : 7.205 
x 


Q ai \ ae 


The above equations may be readily solved on the computer by applying 





the following modifications to Fig. 9, on page 25. 
1. Elimination of Amplifier 5. 
2. Replacement of 4k by 4n when the velocity is negative; 
and by 4n' when the velocity is positive. 
3. Disconnection of dP from Amplifier 1 during the pre- 
soribed interval of combusion, thus making dPp zero. 
Items 2 and 3 above may be accomplished through the use Pete It 
should be noted, however, that since the free piston engine clearance is 
variable, the constant pressure combustion process must be initiated by 
@ change in sign of velocity and terminated when a prescribed amount of 
energy has been added rather than at particular displacements, The 


reason for this is that a change in clearance should not affect the 
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amount of fuel energy added, 

Figure 16 illustrates the computer arrangement to simulate the 
standard diesel cycle. The functions of the additional amplifiers and 
relays are as follows: 

1. Relay 1, actuated by a change in sign of velocity, applies 
the coefficient in when velocity is negative, and the coefficient 4n’ when 
the velocity is positive. 

2. Relay 2, also actuated by a change in sign of velocity, dis- 
connects dP from Amplifier 1 when the velocity goes positive; and connects 
the eirouit 80 that dPp may be reconnected by Relay 3 at the proper time, 


dt 


Se Amplifier 5 performs the integration of 1 n'Pax, whenever 


10 
x is positive. The ee may be rewritten as 
Baer ‘ 7 oF & lfdQ 


This integration, therefore,yields a result proportional to the total 
instantaneous fuel energy. This is compared with the desired value of 
Qe in Amplifier 6 which actuates Relay 3 when {(3) dt equals m¢HV. 
Figure 17 illustrates a typical solution ors this arrangement, It 
should be pointed out that the operation of this analog does not depend 


upon a pre-set stroke or clearance, 
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8. Energy Relationships in the Engine. 

In the operation of the analog of the standard diesel engine, it will 
be necessary to know the net enorgy lost as heat from the system; and the 
net energy added during the combustion process. These quantities may be 
determined by the integration of equations 7.2.2, 7.2.4, and 7.2.6 over 
the appropriate displacements, These are: 


1, Energy as heat transferred during compression 


a= gM 1]. 


2e Net onergy as heat transferred during combustion 


@F a Tore - Vy). one 


Se nergy as heat transferred during expansion 


Qn = arama ao I3 v1 Ft 4)" a 4] mee 


where the subscripts 1,2,5, and 4 refer to the points illustrated in 

Fig. 17. It can be seen that with the computer furnishing a picture 
similar to Fig. 17, all of the required pressures and volumes will be 
known, thus all of these Q's may be determined, The sum of Q and Q: 
will be that portion of Qrej chargeable to the engine; and Qn is an 
idealized mpHV product which does not include heat lost during combustion, 


nor energy lost through incomplete combustion. 
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9, Electronic Analog of a Reciprocating Compressor. 

The force exerted upon the compressor cylinder may be most conven- 
iently represented by referring toa typical indioator diagram for a 
compressor, such as appears in Fig. 4(b), page 9. In the case of the 
free piston device, there is no fixed geometrical point. The only known 
quantities are the pressure limits. 

In a manner similar to that used in the engine equations, the 
differential equation of the first law for the compressor may be written 


as 





ark 


ae fda AQ. _ ; 
a. pe Se eS], 


where A, is the compressor piston area. 
P, is the compressor pressure. 
y is the compressor displacement from the head. 
dQ. is the rate of heat transfer, 


at! 
This equation, however, is valid only where the mass is constant, not 
during the delivery or intake processes, For these processes, it can 
be said that 
Po = Pq, for delivery, and 
9.2 
Po = Py, for intake, 
The problem of the changing mass system of the compressor may be 
readily solved on the computer with essentially the same arrengement as 
that used for the engine, but with the use of diode limiters on the 


amplifier generating Pg. Such an arrangement is shown in Fig. 18, 


which solves the polytropic equation 


die _ Me keg 9.3 
At 
J 
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It oan be seen from Fig, 18 that -P, is always negative. With this 
being the case, the diodes in the feedback of Amplifier 7 will fumction 
as follows: 

l. With y negative Seana ae 8 Soe be positive, and-P,, 
the output of Amplifier 1 , will tend to go more negative. This will 
continue until-P, overrides the bias on Diode 1 at which time this tube 
conduots, thus stopping the integration with-P, held oonste t at-Pa. 
This situation will continue as long as y is negative. 

2. When y changes sign, dP will become negative and the out- 
put of this amplifier will beoome nei negative, When this happens, 
the initial condition will be-Py, whioh has been imposed by Diode l. 

3. The absolute value of P, will continue to deorease until 
the bias on Diode 2 overrides-P,, This tube will then conduct, thus 
holding the output at-Pj. This will continue until y changes sign and 
the integration is started at the initial condition-P;. From this point 
on the cycle will repeat. 

Thus it can be seen that by the use of two diode limiters, equation 
7.1 may be solved in the case of a compressor without taking into con- 
sideration the changing mass problem. Furthermore, this arrangement 
will simulate a reciprocating compressor that has a variable stroke, 
which is precisely the situation in the free piston machine. This is 
because the integrations covering the compression and expansion strokes 
are initiated by a change in sign of the velooity, not at a particular 


position. 


The analog illustrated in Fig. 18 was constructed using the assump- 


tion that the compression and expansion proceed at the same polytropic 


exponent. 
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Figure 19 illustrates two solutions of the compressor problem, 
corresponding to two different stroke lengths. Here it is again noted 


that a pre-set stroke and clearance is not necessary. 
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10. knergy Relationships in the Compressor, 

In order to solve equation 5.5.1 it is necessary to know the amount 
of energy removed as heat from the compressor per cycle, Recall that 
applying the coefficient n,instead of k is equivalent to letting 


a z A. (Fas) 8s ) 10.1 





which produces a reversible polytropic process wherein there is no not 
loss or gain of energy over a cycle. This will be true only in the case 
of a constant mass system. The compressor, however, compresses one mass 
and expands a different mss. By integrating equation 8.1 over the com-~ 
pression and expansion strokes, the net energy transferred as heat per 


cycle is 4 


Qn= - ma RT (K-ne fa) 4] 
Cc K=4)( tea) lA 7 1062 


where mg is the mass of air delivered by the compressor per cycle, This 
equation neglects the energy transferred during the delivery and intake 
processes. It can be seen that if m, is zero, or if k ® no, there will 


be no heat loss. 


*See Appendix III for derivation. 
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ll. Blectronic Analog of a Gas Spring 
The force exerted upon the bounce cylinder is that of a gas spring, 


igs, Pz* = constant. In this case, equation 6.5.1 may be rewritten as 


Gi iin 
dt — z 


where Ph is the pressure on the bounce piston 
zZ is the bounce piston displacement, 


li.l 


On the basis of the preceding analyses, the computer arrangement 
may be immediately depicted as that shown in Fig. 20. The diode in the 
feedback of Amplifier 10 performs the function of an atmospheric relief 
by not permitting Ph to fall below a voltage corresponding to one 
atmosphere, thus simulating the introduction of additional air into 
the bounce soace, A typical solution from this arrangement is illus- 
trated in Fig. 21, which also illustrates the functioning of the ate 
mospheric relief, a. ell 

In the case of the bounce cylinder analog, extreme caution must be 
used in the adjustment of the electronic multiplier forming the Pz 
product, for the reason outlined in section 6.5, This situation is 
aggravated in this case since there my not be an initial condition 
imposed every cycle, That is to say, the solution will be "floating." 
A very slight maladjustment in this multiplier will cause the solution 
to drift rapidly away from the correct value, This malfunctioning is 


illustrated in Fig. 22. This situation must be avoided, 
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Figure 21, Computer Simulation of & Ges Spring. 
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12. Coulomb and Viscous Friction. 
12.1 The Electronic Analog. 

The force of friction is assumed to be composed of the two familiar 
types: Static friction, a constant force always opposing the motion, and 
viscous friction, a force proportional to the velocity, again opposing 
the motion. This total force is illustrated graphically in Fig. 4(d). 
This may also be expressed mathematically as 

m= C-wi (<0) 
Fe ~C-arx ( x >0) y) 


where C is the Coulomb friction force 


TZ. 


if 


v is the viscous friction coefficient, 
The analog to produce a voltage corresponding to this force is 


illustrated below in Fig. 235. 





Figure 23, Electronic Analog of Coulomb and Viscous Friction 

The amplifier has no feedback except for the two diodes and their 
series resistances. Whenever the input is different from zero, even by 
the smallest amount, the amplifier output will tend immediately towards ° 
saturation. This tendency is checked at the voltage "C" by one or the 
other of the diodes, which when conducting, places the feedback resist- 
ance in the cirouit. Figure 24 illustrates the computer analog of the 
force of friction both as a function of velocity and of displacement, 
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12.2 Energy Dissipated by Friction. 
The last remaining item as yet unknown in equation 5.5.1 is that 
portion of Qrej chargeable to friction, 


The force of Friction is given by 


ie = ee 5 > when x is negative, 
and leeeer 
é 
S = eV ; when x is positive, 


In order that this force may be integrated directly, it shall be assumed 


that the free piston has the approximate motion 


Les, Coes 3 cow Ans ot 12.262 


where co, is the engine clearance 
fis the frequency 
Ss is the stroke 


The work done by friction is 


We = § Fy de y 12.2.3 
or om ~ 
. At 
We _ ‘ 5 (SE At. 12.2.4 
O 


Upon substituting and integrating over a cyole, the net energy lost through 


friction work is given by 


We = iG Smee ones 12.2.5 
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13. Computer Operated Relays. 
As outlined in Fig. 16, the use of five relays is required in the 
simulation of a diesel engine, For convenience, these are: 

1, A relay to start and stop integration whenever x is greater 
than x), the port position. This has the effect of putting the "tail" on 
the engine picture, thus simulating tho scavenging process, 

2e Three relays necessary to make dPg zero for a portion of 
the outstroke, thus simulating the constant spereees addition of fuel, 

3e One relay to change exponent from compression to expansion, 
The above relays, of course, must be actuated by quantities available in 
the computer. In all cases the relays wore driven by the output of an 


operational amplifier by use of the scheme illustrated below in Fig. 25. 


w-« ™ 


~ 





Figure 25. Use of Operational Amplifier to Actuate Relays. 
As can be seen above, tho amplifier has no feedback except diodes to 
limit the output to a voltage, €,sufficient to actuate the relay. The 
action of the above arrangement is to yield an immediate step voltage to 
the relay whenever the algebraic sum of the inputs passes through zero. 
Figure 26 illustrates the output of such an amplifier actuated by a 


sine wave, It can be seen that this output is very nearly a square 


WAV o 
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14, Electronic Analog of the Free Piston Generator. 

The preceding sections may be summarized by referring to an overall 
Computer arrangement for the solution of the piston motion. This is 
illustrated in Fig. 27. This diagram is basically the same as Fig. 7, 
except that the details of Figs. 16, 18, 20, 23, and 25 are included, 
The numbers on the amplifiers correspond to those used in the above 


mentioned figures, 
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15. Time and Magnitude Scaling. 
15.1 General, 

All of the preceding sections have been rather qualitative in nature 
in order that the principles of the analog approach to the problem may be 
presented more clearly. In order to obtain a numerical result, however, 
it is necessary to consider magnitude and time scaling within the computer. 
15.2 Scaling the Pressure Analogs. 

As outlined in section 6.5, the various pressure analogs were con= 
structed by using as a driving force a separately generated displaced sine 
wave, operating at a fixed frequency, based on real time. Such will not 
be the case when each of these analogs becomes a part of the overall 
hook-up as illustrated in Fig, 27. In this case, the integrations form- 
ing displacement and velocity will seek some natural frequency, based on 
Computer time, not real time, To this end, a time scaling factor of 100 
will be used. Thus, a real machine operating at 600 cycles per minute 
will appear on the computer as operating at 6 cycles per minute, The 
only criterion necessary in making the pressure integrations compatible 
time-wise with the displacement integrations is that a common time scale 
be used. When this is done, all integrations will praotes at the same 
rate. 

The: magnitude scaling is entirely arbitrary, depending upon the 
maximum value of the variable, and its relationship to the 50 volt limi-~ 
tation of the BEAC amplifiers. The only criterion here is that the same 
scale factor must be applied to the initial conditions, 

There is a definite relationship between the various scale factors 
involved within each computing loop simulating the engine, compressor, 


and bounce pressures. ‘the compressor, for example,solves the equation 


08 








Ge Ze neR& 
R= At |= 16%2.1 


or, in terms of computer voltages 


~~ 
ox ° 
R = a PB OL 15 eZe2 


eadln, AY, 
XL y G 


where CY; is the magnitude scale, units of "1" per volt, and the barred 


GX 


Mode 


quantities represent the computer voltages corresponding to the real 
quantity, By comparing the above equations, it can be seen that Xp, 


may not be arbitrarily assigned, but is given by 


Xp = pe Gan ; 16.2.5 
G OX 


Similar expressions exist in the engine and bounce analogs. These are 


xa = FRX 


F, 15 e204 
é Oy 
and 
OS 2 


As outlined in section 6.2, there is no liberty selecting the 
magnitude scales involved in the multiplications and divisions. This 
scaling has already been performed for all amplifiers except those in- 
tegrating to form the various pressures, 

In the case of the engine, for example, the integrating amplifier 


may be treated as follows: 
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~_ 
Xp 
P. a, ae 
In this amplifier, 
eusts Oo PR. Poe 
= ea tan 
~p, 'E Os d fy? 16.2.6 


where CXy is the time scaling factor, computer time. 


real time 
When the above constants are associated with the circuit parameters as 


outlined in AppendixTT, the following equation is formed 





Gare 
i = oe : 15.2.7 
Xp Ae R, Cs, 


When equation 15.2.3 is substituted, 


a, : Xy 
aa : 15.2.8 
R Cs, Ky 


By similar arguments, the circuit parameters involved in the com- 


pressor and bounce pressure integrations are given bys 





Ay _ Ag 





15.2.9 
Ry Cs fy Xe 
and 
cd X z 
ap C10 Xs Xp , toe @ 
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The subscripts 1, 7, and 10 used in the above equations refer to the 


number of the amplifier as given in Fig. 27. 


Since x, y, and z refer to the samemotion, it must be true that 


Xx = Oy = > > 15e2ei1 
and 
Sy ~ Xo oee te ace 


When this substitution is made, there finally results 


2, a, Aro = Ky 15.2.13 


— 
—— 


R, Csi 7 Ny Csy 7 Kio Cs10 Oy an : 











This relationship specifies the definite time-magnitude relationship 
between the pressure integrations and the displacement integration. 
15.3 Sealing the Displacement Integrations, 

The displacement integrations may best be analysed by redrawing a 


portion of Fig, 27: 


. X 
ae 
=e -X > + x 
al, Qy, 
,* Ib 





x, X 


ff 


= ak Aco, = Acx, P + Ax, R + oF |. 15.301 
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From which 





DZ yy (ae = Ae Cf pe 
Ris M oy : 16.302 
Ana heya _ A & p 
[i oa 
Ang Rows _ Ab Xp, 
FR = no 3; 15.50% 
143 M oly 
and 
Qhrgit Ns 144 = Xe 
— | a aa? 16.3.9 
Rig NM X 
The following may be written for Amplifier 15. 
— Ole a 
xt 
or 
25 OX Xi 
e — ¢ 1oeo6 7 
Ris Csis Oy He 
And for Amplifier 165 
_ NS 
Ci. =e —— dt 
or 
aay oy 
TC = ——xX— , 15.3.9 


By comparing this equation with equation 15.2013, it can be seen that 
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a = 2, _ Ar = _ae = Ai, 
Xx Oe R, Coy yoy eGo a 


This equation means that the RC time constants in the x, Pp, &, and R, 


integrations must be the same, 


15.4 Scaling the Friction Analog. 


Again referring to a portion of Fig. 27: 





15.4.1 
from which 
C13 Keg 2 ICV 
RR - 15.422 
13 nF, 


No further scaling is required in the other amplifiers of Fig. 27, 
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16. The Pescara Machine on the Analog Computer. 
16,1 Results. 

Reference 1 contains the data and a digital analysis of a free piston 
gas generator manufactured by the Pescara organization in France. Appen- 
dix IV contains the applicable original data and some of the results given 
by London. It is interesting to note that the analog approach to the 
problem is nearly the reverse from the digital eppranekt The digital 
approach uses as a starting point an assumed stroke and clearances, and 
then proceeds to calculate some of the geometry, frequency, power output, 
and other quantities of interest. The analog computer approach uses as 
a starting point a particular geometry, and yields as a solution stroke, 
clearance, and frequency, Unfortunately, the data required by the com- 
puter were not entirely complete in the original Pescara paper, 
Conspicuously absent were? 

1. The piston mass, M. 

2e The absolute bounce clearance, b,. 

Se The initial bounce pressure, Phe 
All other items were either given or were determined by simple calcula- 
tion. 

Of course, the computer cannot produce a solution in the absence of 
the above items, London, however, had calculated these items in his 
analysis, and his results being the best available data, were used for 
these three parameters, It should be pointed out, however, that this 
digital analysis is not quite compatible with the computer solution in 
thats 

1, The digital analysis used a modified Otto engine cycle, 
while the computer uses a diesel cycle, 


2 The digital method charges friction individually to the 
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engine and compressor in the form of mechanical efficiencies, while such 


is charged to the machine as a whole in the computer approach. 


Oe 


The digital method uses “lumped” quantities such as net 


engine work, net compressor work, and net friction work. 


Such is not 


possible on the computer since the computer must consider differential 


items in solving the equation of motion. 


With these exceptions noted, the following compares the computer 


solution™ of the Pescara machine with the combined original data and 


London's analysis i 


1. 


Ze 


Geometrical Results 


Item 


Frequency, cycles/min. 


Stroke, inches 

Engine effective stroke, inches 
Engine clearance, inches 
Compressor clearance, inches 
Bounce Clearance 

Mass Results 


Compressor air delivered 
lbs/side/cycle 


Ratio, engine/comp, air 
Aire-fuel ratio 
Gas flow rate, lbs./hr. 


Fuel flow rate, lbs./hr. 


Analog 
600 
17,5 
9.74 
1.10 
Ze08 


254 


0.500 

0.570 
2207 
56,600 


589 


bi 
Refer to Appendix IV for scaling of Pescara Machine 


itResul te and Calculations appear in Appendix V, 
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Pescara 
613 
17.6 
9.70 
1.14 
258 


200 


0.580 
0.516 
5220 

29 ,000 


460 





Se Thermodynamic Results. 


Item Analog Pescara 
Delivery Temperature 1550°R 1405°R 
Alr Horsepower 1200 1138 


The net energy balance wag also calculated, This showed a discrepancy 
of 6.6% between the engine work and the compressor plus frictional work. 
Figures 28, 29, and 30 are the computer solutions of the engine 
pressure, compressor pressure, and bounce pressure, respectively, as 

functions of their respective displacements, 
16.2 Discussion of Results. 

As shown in the preceding section, the geometrical results were 
in excellent agreement with the actual machine data. The rest of the 
results, however, were close, but not as good as might be expected, 
This is attributable to two causes: 

1. The analog computer was actually solving a different kind 
of problem than the Pescara machine, even though the geometry was essene 
tially the same. This is because: 

a, A diesel cycle was used instead of a real engine cycle, 

b. Different thermodynamic constants were used, 

ce. the frictional constants were assumed at much too low 
a value, 

2e Computer error. In this regard it is to be noted that 
extreme difficulty was encountered in making the bounce analog stable, 
This phenomonon resulted from the error in forming the Pz product by an 
electronic device. Since z Changes sign it is imperative that’ 

Pz = ~P(-z), so that the solution will reproduce itself regardless of 


direction of motion. The electronic multiplier is not capable of the 
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~~ | 
1800 ; 
2 tL ate 
HA 


_ S “a 
e) i ie 


Engine Displacement, inches 


Pressure, psia 


' Figure 28, Cemputer Simulation of the Engine of the 
Pesears Mschine, 
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Pressure, psia 


Pressure, psia 


a | 


~~ 
5 oO 
eG 


Ww 
© 
sips 


© 


zs 8 8 


© 


—__ Sae 
2Q 


Compressor Displacement, inches 


Fisure 290 





Computer Simulation cf the Compressor 
of the Pescsre Mzehine, 


Bounee Piston tetaialae a cok 


Ficure 30, 


Computer Sinulation of the Bounce 
of the Pesesrea Machine, 


. 
: 
t 
. ae 
v ’ e 
¢ 





accuracy that is required in the case of the bounce piston. The engine 
and compressor analogs are not sensitive to this effect in that an initial 
condition is imposed each cycle which prohibits the solution from drifting 
AWAY e Thus, the analog of the bounce piston exhibited a phenomenon 
similar to, but not as pronounced as, that showm in Fig. 235, page 51. 
The maximum drift was of the order of 3 volts. The scaling, however 
was such that a 3 volt departure represented a force error of about 
59,000 pounds, which is sufficient to disrupt the entire analog, causing 
the solution to become unstable after about three cycles, Thus the 

| 
solution as outlined in Appendix V approaches, but may not be, the steady 


state, This effect may account for most of the discrepancies, 
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17. Conclusions 

The free piston engine problem may be solved on an electronis 
analog computer, if such a computer is capable of generating a stable 
solution of the gas spring. The computer can produce a solution to the 
equation of the first law of thermodynamics in the simulation of an in- 
ternal combustion engine and of a reciprocating compressor, An analog 
computer is not capable of solving the gas spring problem to the required 
accuracy in this manner due to small inaccuracies in electronic multipli- 
cation. 

Two alternative ways are suggested to produce an accurate analog 
of the gas spring, These ares 

l. Use of a function generator to yield a function which may 
be used to give PzK = constant. Such a device would keep positive control 
of the bounce force at all times, 

2. Use of one-quadrant multiplication. It has been pointed 
out that the primary difficulty was that electronic multiplication could 
not give exactly Pz = -P(-z), This difficulty may be avoided by using 
zZ having always the same sign. This may be accomplished through the use 
of velocity relays and a sign changer to switch both the input and oute 
put of the multiplication such that the multiplier operates only in one 
quadrant, 

It should also be pointed out that the solution of this problem 
required a good sized computer, in that a total of 26 amplifiers, eight 
diode limiters, eight voltage sources, and six eteronic multipliers 
were required, in addition to the necessary read-out equipment, This 
is considered the minimum computer size to handle such a problem. 


Further refinements, or the addition of another compressor analog for 


VP? 
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the solution of the free piston compressor would, of course, inorease 


the computer requiremonts. 
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APPENDIX I 
LIST OF ASSUMPTIONS 
The scavenge air is furnished at constant pressure and 
temperature, 


The delivery to the turbine occurs at constant pressure 
and temperature, 


The compressor intake is 70°F; one atmosphere, 
The operating fluid obeys the ideal gas laws, 


The pressure drop across each valve end engine port is 
five per cent, 


The geometry is symmetrical with respect the centerline, 


The compressor compresses and expands in accordance with 
a polytropic exponent of 1.30. 


The sensible enthalpy of the fuel may be neglected. 


The engine operates in accordance with an air standard 
diesel cycle. 


The energy transferred as heat during the compressor intake 
and delivery is neglected. 


The frictional foree consists of coulomb and viscous 
friction, 


The friction work is determined in advance by assuming 
that the stroke is sinesoidal, 


Other numerical assumptions regarding the specific problem 
of the Pescara machine are found in Appendix IV, 
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APPENDIX MII 


Analog Computer Symbols 


Circuit Symbol 


le Potentiometer: 
&; 


& RK 
Sy = ahh: 
oy 6 ode ie p 
2 2 3 
a3 3 Cy Ce R oz 


€,-oAws 


Se Summing Integrator: 


Cy 
wa 
2 <2. 9G ! 
Qe ? ek SS, 


@, gh i ¢ 


4. Relay Driver: 
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(2) 


Jo BAB, Sar 
Functig 97° - 


CG, = ae, 
¥ 
/ Ra 
+ 29% 
3 é,) 
t 
r 


(a when Lf + 42 >0 
C37 


+b when BE 4 AG O 
K, R2 es 





APPENDIX II (Continued) 


Circuit Symbol Meaning 
5. Beeing Electronic Multiplier: 


Sl Sie 
ah 
4, a [ser >a C3 736 SS 


Note: Direct multiplier output is approximately —- 0, © 5 
5x 10% 


6, Donner Electronie Multiplier: 


: t 
iA e C3= ~ joo 1% 
3 3 


2 


7. Division using Boeing Multiplier: 


& 
«SHO @ 
3 R, e, 





8, Division Using Donner Multiplier: 





Other Special Circuits sare explained in text, 
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APPENDIX III 
Heat Removal from the Compressor 


The compressor cycle is given by the figure below: 


For the processes 1-2 and 3-4, the rate of heat transfer is given by 


da _ A (ie 


p, ZS, 
OG eve oe 


It shall be assumed that the delivery and intake processes involve 


no energy transfer as heat. 


For process l=2 


or 


4a = [A HH) a] 


Upon integrating between the volume limits 1 and A 


= _ ae | a ae 
Gia = Ac (k- (med) ie L sr -. 


- oe Ne~ 
Now Ae = )lne-d) By, pee | - ve 7 
Aey, 2M 5 


PV.= m, RT. 


Therefore 
oe 7 Gear = 
Ora : (-If me - 1) — (T; T) 


76 


and 





By a similar process, Q34 is given by 
. Nain R 
Qs = at, ha, | 73 - 7; 
Meee 2 i 


The net energy transfer is 


OQ. = ), *Gh, = — Ae (1 T)(m,- m)R 


ze (K- Drea la 3 ON 


M37 m, = - Me, 


/ 


Therefore 


Se een Te 
CO). ike pe, ee ee 
a (H- 4) n- 4) A f EE 1 ‘ 


22 (hen) era} ()? = 2| 


(K-1)(ne-4) a 
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APPENDIX IV 
PROGRAMMING THE ANALOG COMPUTER FOR THE PESCARA MACHINE 
XY Data, 
Listed below is the combined data from the original data and the 
results given by Londons (2) 


A, Geometry: 


Ae 0.979 ft (x) 
Ao 5.73 ft° (*) 
Ay 6.77 £6° (x) 
M 2080 1b = (#) 
Xp 10,84" (tT) 
se) 1,24" (tT) 
b 21.14" (a) 
Be Operating Conditions: 
P* 4,85 (*) 
Py 64,.5psia (*) 
PP 293 psia (+) 
C. Assumed Quantities: 
C 50 1b 
v O.l 1b sec/ft 
n 1,30 
= n't 1,40 
7 Ng 1.50 


*Original Data. 
TCalculated from original data. 
+From London's solution. 
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II. Assignment of Soale Pactors. 


A, Pressure Soaling, 


1. 


Ze 


Se 


The Engine, 

Pp = 64,5 psia 

The compression ratio is estimated as 10 
Therefore Peak pressure % 1500 psia 


30 psia/volt = 4320psfe/volt 


Let Xp. 
Then Pp = 2.15 volts. 

The Compressor. 

Py, 3 14.0 psia 

Pq # 14.0 x 4.85 © 67.9 psia 

Let <p = 10 psia/volt = 1440 psfa/volt 


= 1,40 volts 


3) 
” 


- 6.79 volts 


GO 
fi. 
8 


The Bounce. 
Let O¢p = 20 psia/volt = 2880 psfa/volt 


Pe = 14,65 volts 


B. Displacement Scaling. 


Let of, ® 4" per volt = 1/24 ft per volt 


Then! Xp = 21.68 volts 
0 = 2,46 volts 
b = 42.28 volts 


C. Velocity Soaling. 


It is knowmm that 


Xx] + 9 cos 62.8t inches 


Xiax % 50 ft/seo “ 
Let XH = 2 ft/aec/volt 
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+ 
D, Acoeleration Scaling. 

Xnax % 3000 ft/seo® 

let C¢y 2 100 ft/seo*/volt 
E. Time Scaling. 

let ey = 100 seconds computer time/l second real time 
F, Friction Force Scaling. 

let <r # 20 lbs/volt 

then C = 2.6 volts 


G. Summary. 


Sp. = 4320 psfa/volt 
Xp = 1440 psfe/volt 
Sep = 2880 psfa/volt 
Oy = 1/24 ft/volt 
Oly = 2 ft/sec/volt 
Xky = 100 ft/se0"/volt 


Xf, = 20 lb/volt 
Xr = 100 seo/seo 
III, Determination of Circuit Quantities 


A. Pressure and Displacement Integrators (Amplifiers 1, 7, 10, and 
16.) 


From equation 15,5.10: 


2), Qy  _ A | Qiu _ XX 


EE 


RG: Fy Cey ~ RoCso Rye tee 








= 0.480. 


Let R's = 1 meg 


C's & Lat 


ats # 0.480 
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The velocity Integration (Amplifier 15) 


B. 

From equation 15.5.7 

06 ap Se OOo 
Ris Csis Xz Ce 

Let 
Cas = 44S 
ig = 7 meg 
ANG: m0 

C, Summation of Forces (Amplifier 14) 


From equations 15.3.2 to 15.5.5: 


iat Kew = Ae Kee = OG Ge 
Rig) Hfnws 
One Rey Mee 4878 
Riga M oy 
Ary 3 Reig - Ag Op 7 # O22 
Pigs HecGs 
Out Rsie = OK, i [DMO IE. 
Rise Mey; 
Lets 
Reig = 1Lmeg 
nae Eee meg , Thon Gijyy = 0.654 
Risg = +S meq, ten Ayg> 0.637 
(). 30ek 


nis) ig , | meg, Bon A143 
Rivs = 10 meg Then diy = 00310. 
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D, Formation of Frictional Force (Amplifier 13) 


From equation 15.4.2 


21a Ksis = —s . = Or OeL 
RUE Ot r. 

dele 
Rig = O1 meg 


Ry = 10 meq, ther A,3 = 400 
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APPENDIX V 
CAICULATION OF RESULTS FROM THE ANALOG COMPUTER SOLUTION 


Ll. Frequency. 

The measured frequency was six cycles per computer 
minute, coyvwesponding to a real frequency of 600 cycles 
per minute. 
2e sStroke.and clearances. 

Stroke from figure 28 ~- 17.5 inches. 

Engine effective stroke - 9.74 inches. 

Engine clearance - 1.10 inches. 

Compressor clearance - 2,33 inches. 

Bounce clearance - 2.54 inches. 

45 Compressor delivery. 
From equation 5.5.2: mg = Py Ag Yy/ RT, . 
| > - 
Y, = 29.9 volts = 14.95 in. 
P, = 14.0 psia 


Ale 5.73 ft 


To = 530° R 


MG as ate ayy yp z= 0.500 lb 
SeO]( 935 a 


4. Engine charge mass. 


From equation 5.4.5: 
0.7 
m,/Vp * 0.0113 Py 
= 0.209 lb/ft? 


V_2(10.84)(0.979) = 0.884 £t° 
P 12 
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4, (Continued) 


5. Ratio of compressor to engine air. 


Me/ta — 0.185/0.500 = 0.3570. a. 
6. Mass of fuel. 


Frem equation 8.2: mpHV=_k_ Po (V,-V.) 
ns 


B,=51 volts = 1530 psia 
AV = 3.6 volts = 1.80 in. 


= 1.4 (1530)(144)(1.8) = 115,500 ft-1b 
maT = ap aso ace gy 


Assuming a LHV of 18,200 BIU/1b; 


7. Air-fuel ratio. 
aef ratios 0.185/0.00818 = 22.7 
8. Fuel rate. 


Fuel rate=(0.00818) (2) (600)(60) = 589 lb/hr 


Total mass per stroke per side = 0.508 lb. 


Gas rate = (0.508)(2)(600)(60) = 36,600 1b/hr 
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10, Friction work. 


From equation 12.2.5 
)* — 153 ft» 


= (2)(50)(17.5) | ( {od} 
(“re) * (ie 


11. Heat lost by compressor. 


From equation 10.2, (Y. = Thal K- AE ( pr) 7 We- = .* 


(Kr -4)f n-~1) 
Q gz 10.500) 155-5) (530 ip 
1 24) (Oe | (4. @s)° ook -2 | 
= 6170 ftelb. 
l2.Heat lost from engine. 
From equation 8.1, “nes (ken) FV ) EV (8) 
‘ (K-1Xo-1) 
Qe — (0,1)(0.979) (64.5) (144) (10.84 
AG 04) (0.3) (12 (S.5p°e? . 3 7 
8,000 ft lb 


15. Temperature of exit gasses, 
From equation 51 5.5.1: 
Mgh, + mpHV = (me+mg)hy + Ore; 


mga = (0+500)(778) (126.7) 49,200 £t-2b 
mg 115,500 ft~lb 


13,300 ft-lb 


. “er 


“Qre3 
from which 
h, = 151,400 = 382 BIU/1b 
From the tables: 


Ty = 15500R, 
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14. Air Horsepower. 


From the ges tables, an isentropic expansion 


@ 


from the state given in item 13 yields an 


ideal work output of 159 BTU/1lb. 


| AHP .. (139) (778) (0.508) (2) (600) (60) - 1200 
0 5 VOO0 


15. Work balance. 
{ 
{ 
Area under the engine curve = 340 volts”. 
We = (340) (4320) (0.979) (1/24) :=60,000 ft-lbs 
Area under the compressor curve := 160 volts®. 
Wa = (160) (1440) (5.75) (1/24) = 55,200 ft-lbs. 


We = 1535 ft-lbs. 


Discrepancy: 4,650 ft-lbs 
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APPENDIX VI 
THE BOEING SLECTRONIC ANALOG COMPUTSR 

The Boeing Slectronic Analog Computer (BEAC) consists of an equip- 
ment rack, in which is mounted a number of operational amplifiers, 
coefficient potentiometers, limiters, multipliers, and voltage sources. 
The rack also contains a common power supply for the various items, For 
this investigation, two such racks were used, A photograph of this equip- 
ment is shown in Fig. 31. Also illustrated is the x-y plotting device, 
and the time recorder, The hook-up to solve a particular problem may be 
made in two ways 

1. Direotly on the front of the counieee, or 

2. On a pre-patch board, which is then mounted on the front of 
the machine, All the necessary connections are led in beneath the pree= 
patch vanel, 

The Boeing operational amplifier is a high-quality DC amplifier 
having an open loop gain of approximately 50,000. This amplifier will 
operate with stability with any feedback ratio up to unity, It may be 
used for summation, integration, differentiation or more complex 
operations. The amplifier remains linear as long as the output is less 
than 50 Sea, A ape tube oircuit shorts the output at about 51 volts, 
thus this value is an absolute operating limit, 


The general operation of the computer is quite stable, This 
! 
stability bay ke enhanced further through the use of an automatic bale 
ancing system, which restricts the amplifier drift to less than a few 
millivolts per month. 
Detailed descriptions and operating showedures may be found in the 


manufacturer's eta, 
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